arXiv:1502.00242vl [physics.flu-dyn] 1 Feb 2015 


Stokesian swimmers and active particles 

B. U. Felder ho 13 

Institut fiir Theorie der Statistischen Physik 
RWTH Aachen University 
Templergraben 55 
52056 Aachen 
Germany 

(Dated: February 3, 2015) 

Abstract 

The net steady state flow pattern of a distorting sphere is studied in the framework of the 
bilinear theory of swimming at low Reynolds number. It is argued that the starting point of a 
theory of interacting active particles should be based on such a calculation, since any arbitrarily 
chosen steady state flow pattern is not necessarily the result of a swimming motion. Furthermore, 
it is stressed that as a rule the phase of stroke is relevant in hydrodynamic interactions, so that 
the net flow pattern must be used with caution. 
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I. INTRODUCTION 


The dynamics of swarms of active particles has been stndied intensively in recent years 
Q-i- Much of the work is based on the assumption that each particle of the swarm 
moves with a velocity determined by its own activity and by the local fluid flow velocity 
arising from the flow patterns of surrounding particles. The flow pattern of each particle is 
centered on that particle and is carried along with the particle velocity. Due to the perpetual 
change of particle positions this leads to an interesting many-body problem with complicated 
dynamics. 

The assumptions in the theory, as used in practice, can be questioned on two counts. 
First, it is usually assumed that at each point in time the net steady state flow pattern of 
each particle is all that needs to be considered. In fact the net flow pattern must be regarded 
as the time average over a period of the swimming or flying motion. On the fast time scale of 
the period there is an additional oscillating flow pattern. The phase of the oscillating pattern 


is ii^ortant and affects the hydrodynamic interaction and hence the swimming velocities 
[^-0. Second, a net steady state flow pattern is often assumed without derivation from a 
swimming motion on the fast time scale. 

In the following we study the second assumption on the basis of low Reynolds number 
hydrodynamics (^. Therefore the fluid equations of motion are Stokes equations for a viscous 
incompressible fluid, and inertia effects are neglected. In Stokes hydrodynamics the flow at 
each point in space is determined instantaneously by the no-slip boundary condition on the 
surface of each of the particles. 

We study a single distorting sphere and calculate the resulting net flow pattern to second 
order in the amplitude of stroke. It turns out that a commonly assumed active particle flow 
pattern, of so-called B 1 B 2 type, cannot be realized as the result of the swimming motion of a 
distorting sphere. In particular this calls into question the calculation of the hydrodynamic 
interaction of two swimming micro-organisms for which the B 1 B 2 model was first proposed 


We conclude that instead of assuming a particular net steady state flow pattern for an 
active particle it is preferable to consider a swimmer characterized by a combination of low 
order oscillating multipole moments and to calculate the corresponding net flow pattern. 
Several examples of such explicit calculations are presented. Each of the resulting net flow 
patterns can be used in the dynamics of swarms of active particles, though with the caveat 
that the phase of stroke may be relevant in hydrodynamic interactions. 


II. SWIMMING SPHERE 

We consider a sphere of radius a immersed in a viscous incompressible fluid of shear 
viscosity rj. The fluid is of infinite extent and at rest at infinity. It is made to move as 
a result of shape deformations of the sphere, which change the undeformed sphere with 
surface Sq into a body with surface S{t) at time t. The fluid flow equations are formulated 
conveniently in the instantaneous rest frame of the body. It is assumed that in this frame 
the flow velocity v{r, t) and the pressure p(r, t) satisfy the Stokes equations of low Reynolds 
number hydrodynamics 


— Vp = 0, V ■ n = 0. 


( 2 . 1 ) 
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The flow velocity is assumed to satisfy the no-slip condition at the surface S{t). A point 
on the surface Sq of the undeformed sphere is denoted by s, and the corresponding point 
on the surface S(t) is denoted by s -|- with displacement vector The no-slip 


condition reads 11 


v{s + ^{s,t)) = 


dt 


( 2 . 2 ) 


We place the origin of a Cartesian system of coordinates at the center of the sphere Sq. 
By definition 

[ ^dS = 0. (2.3) 


'So 


We also exclude the radial displacement corresponding to uniform expansion of the sphere. 
We assume for simplicity that the displacement is axially symmetric and choose the axis 
of symmetry as z axis. As a consequence the flow velocity and pressure are also axially 
symmetric, and the body acquires a translational velocity U{t) = U{t)ez in the direction of 
the z axis, but no rotational velocity. 

In spherical coordinates (r, 6, ip) the flow velocity can be expanded in terms of a set of 
fundamental solutions of the Stokes equations (2.1), 


with 



CXD 

v{r,t) = -U{t)ez + '^mi{t)ui{r,e) + fc;(t)n;(r, g). 


i=i 


1=2 


/ \ l -\-2 

ui{r,6) = (-] [{I+ l)Pi{cos6)er + Pl{cos6)eg], 

Vi{r,0) = (^^^ [(/+ l)P;(cos^)e^ + ly2p/(cos^)e,)], 


( 2 . 4 ) 


( 2 . 5 ) 


with Legendre polynomials Pi{cos6) and associated Legendre functions P/(cos6') in the 
notation of Edmonds [l^. In the second sum in Eq. (2.4) the term Vi{r,6) is missing on 
account of the requirement that the body exert no net force on the fluid. It follows from the 
expansion Eq. (2.4) that the translational velocity U{t) may be calculated from the identity 

( 2 . 6 ) 

where the integral is over any large sphere centered at the origin and enclosing the body 
completely. The pressure corresponding to Eq. (2.4) is 


p(r, t) =po + 21 ] ^(2/ - l)ki{t) 
1=2 




Pt{cosO), 


( 2 . 7 ) 


where po is the ambient pressure at infinity. 

Provided the sums in Eq. (2.4) converge we can use the expression also for r 
write 


v{r,t)\r=a 


OO CXD 2 

'^Ai{t)Pi{cose)er + '^Bi{t)j^j-^Pl{cos 6)eg, 


a and 
( 2 . 8 ) 
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which defines the coefficients {Ai,Bi] of Lighthill 
(2.8) with Eq. (2.4) we find the relations 

Ai = 2mi — U, Bi = mi + U, 

Ai = (/ + l)mi + (/ + l)ki, Bi = l)"^; 4“ 2 ^^ ~ ^)ki, 

The displacement may be written analogously to Eq. (2.4) as 

00 CO 

^(s, t) = Mi{t)ui{a, 0) + '^ Ki{t)vi{a, 6), 


and Blake [15|]. By comparing Eq. 


(;>i).(2.9) 


( 2 . 10 ) 


1=1 


1=2 


where the term with ni(a, 6 ) is missing on account of Eq. (2.3). 

We note in particular 

ui{r,9) = (-j [2cos6* + sin6* eg] = — (—J + Se^ej. 

Vi{r,6) = - [2C0S6 Cr — Sin6 eg] = -{I + ErBr) ■ Bz, 

3^4 

U 2 {r, 6 ) = [(1 + 3cos26') + 2sin26' 60 ], 

3^2 

V 2 ir,d) = — (1 + 3cos26')e^, 

4r^ 


( 2 . 11 ) 


where I is the unit tensor. The field Ui is identical to an electrostatic dipole field, the field 
Vi is an Oseen monopole flow, the field U 2 is identical to an electrostatic quadrupole field, 
and the field V 2 is a hydrodynamic stresslet or Oseen dipole flow. 

It is convenient to expand the flow velocity v and the pressure p in powers of the dis¬ 


placement ^ as [11 


V = + ..., p = + .... ( 2 . 12 ) 

By expanding the no-slip boundary condition Eq. (2.2) we find that the velocity at the 
undisplaced surface is given by 


LLg — ^ \r=a — 


dt' 


= V^‘^^\r=a = 


The translational velocity U{t) has the corresponding expansion 

u = + .... 


(2.13) 


(2.14) 


Here the first order term is missing on account of Eq. (2.3). From Eq. (2.9) this implies 
= 2m^i'^ and B^^ = The second order term in the velocity is given by 

1 


= 


47ra^ 


'So 


■u® dS. 


(2.15) 


For periodic displacements with period T = 271/u we put 

Miit) = a{pis cosut — file sin ut), 
Ki (t) = a{hiis cos oot — nic sin ut ), 


(2.16) 
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with dimensionless coefficients fiis, fiic, i^is, i^ic- Then we have from Eq. (2.13) for the hrst 
order flow velocity 


V 


( 1 ) 


(r, t) = —au 


1=1 1=2 


(2.17) 


with multipole coefficients 


= file cos ojt + hzsSincut, 

Ki{t) = Kic cosut + Kis sinut, (2-18) 


The velocity and the rate of dissipation 


in the notation of Felderhof and Jones 

to second order in the displacement can be expressed as bilinear expressions 
terms of the coefficients {fiis, i^is, ««c}- 

From Eqs. (2.4) and (2.17) we hnd for the hrst order moments 


12 


m 


In particular, from Eq. (2.9) 

= 


-aufii{t), 

= —auKi{t). 

(2.19) 

-2aufii(t), 

= —aufii{t), 

(2.20) 


since = 0. 


III. NET FLOW PATTERN 

In this section we consider a periodic swimmer with hrst order how velocity given by Eq. 
(2.17), and calculate the mean second order how pattern. The mean is calculated as the 
time average over a single period T = 27r/u. Thus we consider 

d{2)(j.) = ^ f dt (3.1) 

T Jo 

corresponding to some stroke or displacement Here the how pattern is 

dehned in the instantaneous rest frame at time t and is the solution of the Stokes equations 
(2.1) which tends to at inhnity and has boundary value at r = a given by Eq. 

(2.13). 

The second order time-averaged how pattern may be expanded as 


u(2)(r) = — -I- ^ mp^'U/(r, 6*) -I- ^ A:p^u;(r, 6*), (3.2) 

1=1 1=2 

corresponding to Eq. (2.4). At the surface the mean second order how pattern is given by 

vi‘^) (r) (s) = • V'u(b \r=a- (3.3) 

The right hand side can be calculated for a given stroke ^{s,t). The how pattern v^'^'>{r) 
in Eq. (3.2) tends to —at inhnity with value given by Eq. (2.15). We denote the 
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corresponding coefficients given by the analogue of Eq. (2.8) as {A[^B'j} = {Af'\ . 

These may be calculated from u^g\s) by using the orthonormality relations of the Legendre 
functions 13|. It may be checked that A'^ and B[ satisfy the relation 


Ui^) = -{2B[-A[). 

We dehne the corresponding net flow pattern as 


(3.4) 


v'{r) = t;(2)(r-) + f/( 2 )e 2 


(3.5) 


This tends to zero at inhnity at least as fast as 1/r^ and can be identihed with the flow 
pattern of an active particle. Conversely the question arises as to whether a chosen steady 
state flow pattern can be identihed with the net how v'{r) of a periodic swimmer. For 
example, can we hnd a stroke ^{s,t) for which all coefficients {A'i\ vanish and only B[, B 2 
diher from zero? 


IV. SIMPLE SWIMMERS AS ACTIVE PARTICLES 

In this section we consider some simple swimmers with strokes characterized by combi¬ 
nations of low order multipole moments. The analysis suggests the how pattern of corre¬ 
sponding active particles. 

The calculation of the time-average u^g\s) in Eq. (3.3) is performed most easily by using 
complex notation 


Hi Hic T ) 

nf = Kic + iKis, Ki{t) = (4.1) 


with the identity 

with complex 
^ given by 


u 


( 2 ), 


1 




(4.2) 


V C V — / 

s V / 2 

given by Eq. (2.17) with complex coefficients {niit), Ki{t)}, and complex 


^ = —ia 


1=2 


1=1 


—iojt 


(4.3) 


We have checked the expressions for the i?(}-coefficients given below by a separate 
calculation of the coefficients {A\‘^\t), Bl‘^\t)} and a subsequent time-average. 

We consider hrst a simple swimmer with only Hit ^ 2 ^ H 2 diherent from zero, corresponding 
to the superposition of a potential dipole held, an Oseen dipole, and a potential quadrupole 
held. The mean swimming velocity is derived from Eq. (2.15) as 


7/(2) _ aUj{HlcK2s — hls^2c + ^HIcH2s ~ ^HlsH2c)-i 

5 


(4.4) 
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in agreement with Eq. (7.8) of Ref. 12. From Eq. (3.5) we find for the A', R'-coefficients 


(1Uj(^^IqK2s /^1s^2c “1“ 3^1c^2s 

5 


3/ils/^2c 


B; = 2^(2) + -A[, 

2 9 

^2 = ^^2 = - y aUj{lJ,2cK.2s 


h2s^2c), 


= 10 f/( 2 ) _ 1624), R' = 21 f/( 2 ) _ 3324(, 

a: = ?b; = 6A' 


(4.5) 


The coefficients for I > 4 vanish. The corresponding multipole moments {m(, fc;} are 
calculated from the inverse of Eq. (2.9). The lowest order moments are 


m 


; = + a;), 


U' — _R' 

1^2 — g 2^2) 


rrin = -BL 
2 3 2 


(4.6) 


Squirming at / = 2 with implies A '2 = 0, B '2 = Q and . In the notation 

of Drescher et ah [1^ the source doublet strength is Agd = 3m[a^ and the stresslet strength 
is Astr = 1 ^ 20 ^. These authors measured the values for swimming Volvox, but they did not 
hnd a contribution from the quadrupole m 2 , or higher order multipoles. 

For the mean rate of dissipation or power we hnd from Eq. (7.15) of Ref. 12 


T>(2) _ + /r^^) + 9{kI^ + kI^) + + p^^) + 24(K2ch2c + «(2sh2s)] • (4.7) 

The calculations in Eqs. (4.4) and (4.7) are performed most easily by using the expressions 
given in Eqs. (7.11) and (7.17) of Ref. 12. From Eqs. (4.4) and (4.7) we derive an expression 
for the swimming efficiency dehned by 


Et = 4:r]ua^ 


Wo 


(4.8) 


Without loss of generality we can choose the phase such that pis = 0. We then hnd that 
the efficiency is maximal for 


^2c — 0) h2c — 0, 


K-2s = “ 3 ^ hlc 


h2s — 


11 

5\/2 


hlc 


(4.9) 


with value Et = 5 /( 67 r\/ 2 ) = 0.188. The net how pattern v'(r) can be calculated from the 
multipole coefficients {m(, /c(} by using Eqs. (3.2) and (3.5). We note that the particular 
coefficients ^42 B 2 , as well as and R 4 , vanish when Eq. (4.9) holds. It follows 

from Eq. (2.9) that then k '2 and m'^, k'^ also vanish. The net how corresponds to 
a potential dipole of strength m\ = ( 59 \/ 2 / 75 )/i 4 ^aa; and equal multipoles at / = 3 with 


fcg = m'g = ( 459 / 472 ) 7714 . The measurements of Drescher et ah [1^ for Volvox do not 


correspond to the pattern for optimal swimming in the above sense, since they hnd Agt^ to 
be diherent from zero. 

In order to visualize the axisymmetric how pattern it is useful to introduce a Stokes 
stream function 7 / via the relations T3 


Vr = 


1 d'4> 

r2 sin 9 d9 ’ 


ve = 


—1 dip 
r sin 6 dr 


(4.10) 
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The uniform flow e.z corresponds to the stream function 


^ ^^sin^^, (4.11) 

and the flow patterns in Eq. (2.5) correspond to 

iJui{r,0) = sin 6'P/(cos 6^), ^/’^,^(r, 6*) = sin (cos 6^). (4.12) 

The streamlines of the flow are given by lines of constant ip. In Fig. 1 we show the 
streamlines of the net flow v'{r) calculated from and the optimal moments of Eq. (4.9). 
In Fig. 2 we show the values of and Vq/U^‘^'> at r = a as functions of the polar angle 

6. In Fig. 3 we show the values of ^ at r = a, r = 1.25a, and r = 1.5a as functions 

of the polar angle 6. 

As a second example we consider a squirming swimmer characterized by coefficients 

= 0, ^2 = —/i 2 ) ^3 = “hs with all higher order moments vanishing. The mean swimming 
velocity is derived from Eq. (2.15), or from the expression in Eq. (7.11) of Ref. 12, as 


^ aa;(/i2c/^3s - h2syU3c)- (4.13) 

35 

For the mean rate of dissipation we hnd 

p(2) = 8TT7]u^a^^ + /i^) + 8{f4c + /^D] • (4-14) 

The calculation based on the analogue of Eq. (2.8) for the boundary value u^\s) of the 
flow pattern v(^^(r) shows that all A'-coefficients vanish and it yields for the P'-coefficients 

b; = I m, = 0, 

7 _ 25_ 

= s; = 0, b; = --C/<2), (4.15) 

The corresponding multipole moments are found from Eq. (2.9) as 

m\ = ^ P(2)^ = A;' = 0, 

"^3 = “^3 = ^ m '^ = k '^ = 0 , 

‘2^ _ 

"^5 = -K = ( 4 - 16 ) 

The moments for / > 5 vanish. The net flow pattern v'{r) can be calculated from the 
multipole moments by using Eqs. (3.2) and (3.5). 

Without loss of generality we can choose the ph ase s uch that /r 2 s = 0. We then hnd that 
the efficiency Et is maximal for /i 3 c = 0, /i 3 s = ^/9/8i.i2c with value Et = 12\/2/(357r) = 
0.154. In Fig. 4 we show the streamlines of the net how v\r) calculated from and 



the set of coefficients {A[,B'i} given by Eq. (4.15) for the optimal moments. The net flow 
pattern at the surface r = a is given by 


v'\r=a = 


35, 


cos 9 Br + r — sin 6 ^ H—-(7 + 5 cos 29) sin^ 6*1 eg 

L 54 j 


f/( 2 ). 


(4.17) 


The radial component arises from the second term in the hrst line. In Fig. 5 we show the 
values of v'^jU'A) and v'^jU^A) at r = a a.s functions of the polar angle 9. In Fig. 6 we show 
the values of gX r = a, r = 1.25a, and r = 1.5a as functions of the polar angle 9. 

It follows from Eq. (4.15) that the squirmer can be identihed with a particle 

with particular ratios of the coefficients. 

As a third example we consider a swimmer characterized by coefficients K 2 , with 
all other moments vanishing. The mean swimming velocity is derived from Eq. (2.15), or 
from the expression in Eq. (7.11) of Ref. 12, as 


= — au{7nicf^2s - 7filsl^2c + GK2 cK3s “ SKss^c)- (4.18) 

For the mean rate of dissipation we hnd 

W = [630(pL + hL) + 567{kI + 4) + 1180(4 + 4)] • (4-19) 


The calculation based on the analogue of Eq. (2.8) for the boundary value u^p{s) of the 
flow pattern v^‘^\r) yields for the non-vanishing A'- and R'-coefficients 

A[ = ij?; = m, 

, 2 , 20 

A 2 5 ^^ ®^(hlc^3s hls^Sc)) 

1,6 

A 3 2 ^^ 5 3/iic/t2s + 3p.i5/t2c T ^2c^3s ^2s^3c)) 

2 

A' — _R' — -A' 

/I 4 — ~ ^2i 

1 30 

A's = 2^5 = y - K2sK3c)- (4.20) 


Without loss of generality we can choose the phase such that = 0. We then hnd that 
the efficiency is maximal for 


K 2 c = 0 , 


5 /230 


'''*“3V413 


27 

^3c yTT+lc) 

59 


= 0 . 


(4.21) 


with value = ■\/ll5/826/(37r) = 0.040, showing that this type of swimming is rather 
less effective than for the hrst two examples. For the optimal swimmer the coefficients 
A 2 , B 2 , A 4 , R 4 vanish. 

have considered a similar swimmer with added potential hows 
In our analysis [lH the {mi}- and {fc^j-multipoles can 


Ghose and Adhikari 

3 LC 

r /bo 


18 


with = 


and m 3 = —f/c 3 
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have independent values. We recall that the above expressions are calculated from integrals 


bilinear in the hrst order flow velocity. The work of Ghose and Adhikari [1^ suggests that 


the time-dependent swimming velocity and flow pattern must be calculated to third order 
in the displacement ^ in order to agree with experimental observations of the swimming of 
Chlamydomonas [l^. Delmotte et ah have studied the same swimmer in computer 
simulation. 

We note that it follows from Eq. (2.7) that in all three examples the net flow pattern 
is accompanied by a steady pressure pattern p'{r). In the hrst and third examples the 
coefficients A 2 and B 2 are related by A '2 = \B 2 , and in the second example both coefficients 
vanish. It is not possible to design a swimmer with surface displacement (4.3) for which the 
steady state net second order how pattern v'[r) has all coefficients {A(} vanishing and only 
B[,B 2 of the {i?(}-coefficients nonvanishing. 


V. DISCUSSION 

The net steady state how pattern of a periodically distorting sphere can be calculated for 
a chosen set of oscillating multipolar how patterns determined by the stroke. As we have 
shown, such a calculation yields a net how pattern characterized by a set of steady state 
multipole moments, as exemplihed in Eqs. (4.6) and (4.16). In the study of hydrodynamic 
interactions between swimmers it is preferable to start from a particular set of oscillating 
multipole moments, rather than from a chosen steady state how pattern. The net how 
pattern must be used with caution since the hydrodynamic interactions may be ahected by 
the relative phase of swimming strokes. 

The explicit calculations of Sec. IV provide examples of simple multipolar swimmers. 
We suggest that the corresponding how patterns be used in the study of the dynamics of 
swarms of active particles. It may be preferable to use the stroke of optimum swimming 
efficiency within the chosen class of strokes. 


10 




[1] S. Ramaswamy, Ann. Rev. Cond. Matt. Phys. 1, 323 (2010). 

[2] T. Vicsek and A. Zafeiris, Phys. Rep. 517, 71 (2012). 

[3] M. C. Marchetti, J. F. Joanny, S. Ramaswamy, T. B. Liverpool, J. Prost, M. Rao, and R. A. 
Simha, Rev. Mod. Phys. 85, 1143 (2013). 

[4] D. Saintillan and M. J. Shelley, C. R. Physique 14, 497 (2013). 

[5] J. Elgeti, R. G. Winkler, and G. Gompper, arXiv: 1412.2692. 

[6] G. M. Pooley, G. P. Alexander, and J. M. Yeomans, Phys. Rev. Lett. 99, 228103 (2007). 

[7] D. Saintillan, Physics 3, 84 (2010). 

[8] B. U. Felderhof, Eur. J. Mech. B 31, 168 (2012). 

[9] J. Happel and H. Brenner, Low Reynolds number hydrodynamics (Noordhoff, Leyden, 1973). 

[10] T. Ishikawa, M. P. Simmonds, and T. J. Pedley, J. Fluid Mech. 568, 119 (2006). 

[11] B. U. Felderhof and R. B. Jones, Physica A 202, 94 (1994). 

[12] B. U. Felderhof and R. B. Jones, Phys. Rev. E 90, 023008 (2014). 

[13] A. R. Edmonds, Angular Momentum in Quantum Mechanics (Princeton University Press, 
Princeton, N.J., 1974). 

[14] M. J. Lighthill, Gommun. Pure Appl. Math. 5, 109 (1952). 

[15] J. R. Blake, J. Fluid Mech. 46, 199 (1976). 

[16] K. Drescher, R. E. Goldstein, N. Michel, M. Polin, and 1. Tuval, Phys. Rev. Lett. 105, 168101 
( 2010 ). 

[17] D. J. Acheson, Elementary Fluid Dynamics (Clarendon, Oxford, 1990). 

[18] S. Ghose and R. Adhikari, Phys. Rev. Lett. 112, 118102 (2014). 

[19] J. S. Guasto, K. A. Johnson, and J. P. Gollub, Phys. Rev. Lett. 105, 168102 (2010). 

[20] B. Delmotte, E. E. Keaveny, F. Plouraboue, and E. Climent, arXiv:1501.02912. 


11 



Figure captions 


Fig. 1 

Streamlines of the net flow v'{r) for /ii^ = 1, fiu = 0, other moments given by Eq. (4.9), 
and vanishing higher order moments. 


Fig. 2 


Plot of the components (drawn curve) and Vg/U^‘^^ (dashed curve) of the net flow 

at r = a as functions of 6 for moments corresponding to Fig. 1. 


Fig. 3 


Plot of at r = a (drawn curve), r = 1.25a (long dashes), and r = 1.5a (short 

dashes) as functions of the polar angle 6 for moments corresponding to Fig. 1. 


Fig. 4 

Streamlines of the net flow v'{r) for file = 0, fiu = 0, /i 2 c = —^ 2 c = 1, h 2 s = i^ 2 s = 0, 
fJ'Se = = 0, fiss = ~f^ 3 s = a/9/8, and vanishing higher order moments. 


Fig. 5 


Plot of the components (drawn curve) and Vg/U^'^'^ (dashed curve) of the net flow 

at r = a as functions of 6 for moments corresponding to Fig. 4. 


Fig. 6 


Plot of ^ at r = a (drawn curve), r = 1.25a (long dashes), and r = 1.5a (short 

dashes) as functions of the polar angle 6 for moments corresponding to Fig. 4. 
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